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We analyze the electromagnetic amplitude for the leptonic decays of pseudoscalar 



> 

\f^ . mesons in the quark model, with particular emphasis on r; ^ /^/^ {I = e, ^). We eval- 

^^ I uate the electromagnetic box diagram for a quark-antiquark pair with an arbitrary 

distribution of relative three-momentum p: the amplitude is obtained to all orders 
CN ' in p/ruq, where niq is the quark mass. We compute Bp = T(rj — > l~^l~)/T(rj -^ 77) 

using a harmonic oscillator wave function that is widely used in nonrelativistic (NR) 



quark model calculations, and with a relativistic momentum space wave function 
that we derive from the MIT bag model. We also compare with a quark model calcu- 



Q I lation in the limit of extreme NR binding due to Bergstrom. Numerical calculations 

i-^ ' of Bp using these three parameterizations of the wave function agree to within a few 

percent over a wide kinematical range. Our results show that the quark model leads 
^ . in a natural way to a negligible value for the ratio of dispersive to absorptive parts 

of the electromagnetic amplitude for ry — > /U^;U~ (unitary bound). However we find 
substantial deviations from the unitary bound in other kinematical regions, such as 
r],TT^ -^ e"'"e~. Using the experimental branching ratio for 7/ — > 77 as input, these 
quark models yield B(ri -^ fi^fi^) ~ 4.3 x 10^^, within errors of the recent SAT- 
URNE measurement of 5.1 ib 0.8 x 10^^, and B{r] -^ e^e^) ^ 6.3 x 10^^. While an 
application of constituent quark models to the pion should be viewed with particular 
caution, the branching ratio B{7r^ — > e^e~) ~ 1.0 x 10~ is independent of the details 
of the above quark model wave functions to within a few percent. 



I. INTRODUCTION 

The rare leptonic decays of pseudoscalar mesons, such as tt" -^ e+e^ and ri,K ^ l^l~ 
{I = e,fi), provide sensitive probes of new physics both in and beyond the standard model 
lU . In the case of the vr" and r] decays, a careful analysis of the electromagnetic contribution 
to the amplitude is required in order to isolate a possible contribution from new physics. 
While the absorptive part of the electromagnetic amplitude can be reliably determined by 
unitarity from experimental data on the two-photon width of the pseudoscalar, theoretical 
estimates of the dispersive part are in some disagreement . 

An improved experimental measurement of the branching ratio for rj — > yU+yU^ has 
recently been obtained at the rj meson facility SATURNE in Saclay [^, with the result 
B[rj — >■ /i'''/i~) = 5.1 ± 0.8 X 10^®. This is to be compared with the lower limit obtained 
by neglecting the dispersive part of the electromagnetic amplitude, yielding the so-called 
unitary bound B{ri -^ fi^f^~) > 5^"'* = 4.3 x lO'^. 

The SATURNE measurement is significantly closer to the unitary limit than previous 
experiments (for a compilation of earlier measurements see Ref. |Q). It is therefore of 
interest to reconsider the theoretical situation with respect to the magnitude of the dispersive 
contribution to the electromagnetic amplitude in this general class of decays. 

The unitary bound for the leptonic decay of a pseudoscalar meson V is most conveniently 
expressed in terms of the ratio of leptonic to two-photon widths 
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where mp is the meson mass, and v is the lepton velocity in the center-of-mass: 
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To establish our notation, we express R in terms of a ratio of the amplitudes for the two 
modes. The invariant amplitude for the leptonic decay can be parameterized as 



e^ 



MiV^ l^n = -—mi uCk)Yv{k)L, (3) 

where the explicit factor of the lepton mass mii reflects the helicity flip that occurs over the 
lepton line in Fig. |I[ The width is given by 
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T{V -^ l+r) = —m]mpv\L\'^. (4) 

The amplitude for the two-photon decay takes the form 

M{V^ 77) ^ teh,,^pe^,e'',q1q^F, (5) 

where £1^2 and gi 2 are the polarizations and momenta of the two photons. The form factor 
F is purely real for the decay to on-shell photons. The width is 



«2 



T{V ^ 77) = -j^mlF\ (6) 

and the form factor R for the relative branching ratio Bp of Eq. (|T|) follows 

Unitarity implies a connection between ImL and F, resulting in the following model- 
independent result for the absorptive part of i? p 

Imi?=^lnfl^V (8) 

The unitary bound on Bp is obtained by assuming that the dispersive part of R is negligible. 
In the case of r/ — >■ n~^n~ , Eq. (|^) implies Bp > B^^^ = 1.1 x 10~^. When combined with the 
experimental result for the two-photon branching ratio (Ref. [Q), this leads to the unitary 
limit for B{r] —>■ ^'^ ^~) quoted above. In general, a nonvanishing dispersive part leads to a 
correction to the unitary bound given by 

B- -,^(^\\ (9) 
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Most theoretical calculations of the dispersive part of R have been based on a point- 
like coupling of the pseudoscalar to off-shell photons , including dispersion relations and 
vector meson dominance. Alternatively, a point-like coupling of the pseudoscalar to nucleons 
or quarks (which then decay to virtual photons through a triangle diagram) has also been 
considered |]^. Recently, a calculation of r] and vr meson leptonic decays has been made in 
chiral perturbation theory [Q. 

We consider instead a bound state description of the quarks which comprise the meson. 
We evaluate the leptonic decay V — > l^l^ at the quark level, which proceeds through the 
one-loop diagrams illustrated in Fig. 0. This approach was considered by Bergstrom in 
the limit of extreme nonrelativistic binding, where the bound state quarks are assumed to 
be at rest in the center of mass of the meson (and the meson mass is assumed to be exactly 
twice the quark rest mass). These approximations must be viewed with caution however 
when applied to light mesons such as the rj. 

In this paper we allow the quarks to have an arbitrary distribution of relative three- 
momentum p. We compute the quark model amplitude for V -^ l^l~ to all orders in p/rrig, 
where mq is the quark mass. A quark model wave function supplies the distribution of 
relative three-momenta for the quark-antiquark (qq) bound state. 

This approach is well known in nonrelativistic (NR) quark models ||T0|-|T3|, and has been 
used in a variety of applications, including calculations of the two-photon widths of light 
pseudoscalars |jl^,|12|. We make use of two parameterizations of the momentum space wave 



function: a harmonic oscillator form that is widely used in NR quark model calculations 



12,13], and a new relativistic parameterization that we derive from the MIT bag model. We 
also compare with the extreme nonrelativistic limit obtained by Bergstrom. 

Although the NR quark model and the MIT bag model both give successful phenomeno- 
logical descriptions of light hadrons, we use these models here mainly because they provide 



simple analytical expressions for the momentum space wave function in two limits that 
characterize a wide class of models. Our purpose in this paper is to make an estimate of 
momentum-dependent effects in the quark model amplitude for V — >■ /^/~, rather than to 
assess the the predictive power of a particular model. 

We expect that some of the model dependence inherent in a description of light quark 
binding may cancel in the relative branching fraction Bp. In fact, our numerical calculations 
using the above parameterizations of the wave function agree to within a few percent over 
a wide kinematical range. Our results show that the quark model leads in a natural way 
to a negligible value for the ratio of dispersive to absorptive parts of the electromagnetic 
amplitude for rj -^ fi^fi^. On the other hand, we find substantial deviations from the unitary 
bound in other kinematical regions, such as rj, 7r° — > e^e~ . 

The rest of this paper is organized as follows. In Sec. II, we describe our method for 
the evaluation of P — > l^l~ for a quark-antiquark pair with an arbitrary distribution of 
relative three-momentum. We evaluate the electromagnetic box diagram for the leptonic 
decay in closed form. In Sec. Ill, we present quantitative results using a harmonic oscillator 
wave function. In Sec. IV we derive a relativistic momentum space wave function based on 
the MIT bag model, which we also use to obtain quantitative results. We summarize our 
findings in Sec. V. 

II. "MOCK MESON" METHOD 

A. General framework 

A conventional approach to the evaluation of hadronic matrix elements is to decom- 
pose the bound state into a superposition of free plane wave quark-antiquark (gg) pairs 
I0|-p!2[ . An economical and physically reasonable description of the momentum space wave 



function is obtained by assuming that the quark and antiquark have equal and opposite 
three-momentum in each component of the plane wave expansion. In this "mock meson" 
description, the momentum-space state vector \M[P)) for a meson in the center-of-mass 
frame [P = {rrip, 0)] has the following decomposition |T^ 

|M(P)) = v/2^|^^<|.(p)^|g(p))|g(-p)), (10) 



where Ep = wp^ + m^, and rriq is the quark mass. We have omitted color, spin, and flavor 
indices in Eq. (p^ ) for convenience. We use invariant normalizations throughout this paper, 
e-g-, {q{p')\(1{p)) = '2Ep(27i)^S^{p' — p). The wave packet amplitude $(p) is normalized 
according to 

j£pMp)\' = i. (11) 

For the ground state pseudoscalar mesons of interest here, we assume that the wave function 
is spherically symmetric, $(p) = $(p). 

The right hand side of Eq. (0) is a (zero) momentum eigenstate, by construction. 
However, since the quarks are taken to be on-shell, the energies of the individual qq plane 



wave components are in general not equal to the bound state energy. This energy "mismatch" 
leads to several ambiguities in mock meson calculations, including the value to be used for 
the total qq energy running through intermediate states in the amplitude, and in phase 
space factors ||l2|,|l3l. A popular prescription is to identify the total qq energy appearing in 
the amplitude with the mean energy of the wave packet ||T2[. In the case of quark model 



calculations of the pseudoscalar two-photon width for example |jT^ , this prescription leads 
to a phase space dependence on the meson mass that is in agreement with phenomenological 
estimates of the r] — rj' mixing angle (see, e.g., Ref. 0]). 

For the purpose of calculating the amplitude R however, it is crucial to take the total 
energy of the wave packet running through intermediate states to be equal to the physical 
meson mass (the relative branching fraction Bp is very insensitive to overall phase space 
factors, on the other hand). This is necessary in order to obtain the correct unitarity relation 
Eq. (I) between the absorptive part of the amplitude for the leptonic decay, and the (on- 
shell) two-photon matrix element. This prescription has also been used in other "mock 
calculations 



meson 
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In the extreme nonrelativistic limit Eq. (|T0|) implies that the qq annihilation amplitude is 
proportional to the coordinate-space wave function (or its derivatives) at the origin |lTO| , p!T| . 
This limit has been used for example to make calculations of heavy quarkonium matrix 
elements |jl6[|, and was used by Bergstrom in a calculation of P — > /+/~ 0. 

Equation (|10D has been used to all orders in the qq relative momentum in nonrelativistic 
quark model calculations of various matrix elements [|12|,|13] , including pseudoscalar meson 
decay to two photons, Eq. (^). The form factor F is found by evaluation of the Feynman 
diagrams in Fig. | p|,|12[ 



4^3 



■m„ 



271 nip 



Q^ / pdp^{p) In 



Ep + p 
Ep-p^ 



(12) 



where the logarithm comes from an integration of the quark propagator over angles (assum- 
ing that the wave function $ is spherically symmetric). Q is the quark charge in units of 
the proton charge. For decays of flavor-mixed states such as the r], it is understood that 
an expectation value of the above expression for F is taken between the meson flavor wave 
function and the vacuum. 



Equation (12) is obtained by assuming that each plane wave component of the qq wave 
packet has the same total energy |]T1|,|T2[. In our case, this means that a factor of 1/mp is 
extracted from the quark propagator, which is contained in the overall factor of the meson 
mass in Eq. (O). 



B. Application to leptonic decays 



We derive the mock meson amplitude for V — > l^l along the same lines which lead to 
Eq. ([I2D for the two-photon matrix element: 



M{V ^ i+r) = Vh^ I 



d^p 



(27r)3/2 2Ep 



Hp)M,g{p) 



(13) 



where Aiqq is the amplitude for a given three-momentum component of the qq wave packet 



M,,{p) = -i2V?>e'Q^ j -^^Q^^{q)L,M 



{27ry 

[q"^ + ie][{P - qy + ie][{p - qy - m^^ + ie][{q - ky - mf + le]' 

The factor of 2 above accounts for the equal contribution of the two diagrams in Fig. ^ and 
a/3 is a color factor. Q^'^ is the spin-singlet qq current 



1 4? 

Q>^^^ — Y: Xv-xi-ph" [^- ^+m,] ^^uxip) = -^m,e°"^^g„, (15) 



where, e.g., U\{p) is a positive energy spinor of three- momentum p and angular momentum 
A/2 along a fixed quantization axis. L^^ is the spin-singlet projection of the lepton current 

Lt^- = 2i^e°"'^'^g„ u(kh\(k) (16) 

mp 

[spin-triplet components of L^'^ vanish under contraction with Q'^'^ , or under integration in 
Eq. dg)]. 

Comparison of Eqs. ([T3|)-(p!6D with Eq. (|]) yields the following expression^ for the form 
factor L 

Ulq 2 f '^'^P - ' ^ 1 



where 



^-^«^7ffe«V(2^*<''>^^(P'- ('^' 



^^•^^ ^m^J"^ ^ [g2 + le] [{P - qf + ie] \{p - qf - m\^ %e\ [(g - kf - mf + le] ' ^^^^ 

This integral can be evaluated analytically. An identity relating q^ to a linear combination of 
(inverse) propagators reduces the integral to a sum of scalar three- and four-point functions, 
plus an integral of three propagators with a factor of q^ in the numerator. A fictitious 
photon mass is introduced in intermediate calculations as the scalar vertex and box functions 
obtained in this way are infrared divergent (the total integral is infrared finite). We evaluate 
the divergent integrals using expressions provided in Ref. |]T^. After some algebra, we find: 



Ap) = . . ,.! . = \Hxp)lL + Ir] , (19) 



4w (j9 ■ ky — m^m 



^mf 



where 



^ +21„(l-4)-- 



lL = 2 + \n{ -^ + 21n(l - xl) - - \n{xp) + in, (20) 



^As with Eq. (^) for the on-shell two-photon form factor, it is understood that an expectation 
value of Eq. ( |l7|) for L is taken in the case of a flavor-mixed state. 




and 



+ Sp(a;5) + Sp 1-a;,^ +Sp(l-a;,^), (21) 



p ■ k — J {p ■ ky — rn^mf 



Sp(x) is the Spence function 



Sp(x) = - r dt ^""^^ ^\ (23) 

JO t 



A crucial intermediate step in our evaluation of the loop integral in Eq. (p!8| ) is the 
identification of the total energy of each component of the qq wave packet with the physical 
meson mass (i.e., 2Ep = uip). This is necessary in order to satisfy the unitarity relation 
of Eq. (^. In particular, without this prescription for the wave packet energy, the loop 
integral would acquire unphysical branch cuts that do not correspond to the "unitarity cut" 
through the intermediate photons in Fig. 0. On the other hand, we do use the actual plane 
wave energy p^ = Ep = Jp"^ + m^ in our final expression Eq. (0) for /(p). We note that a 
similar prescription has been used in other mock meson calculations, such as the two-photon 



width (cf. Eq. (^ and Refs. HJl 



The fact that our final expression for Im L exactly satisfies the unitarity relation of Eq. 
(§) is a nontrivial check of the above prescription for handling the ambiguity in the wave 
packet energy. Indeed, after an analytical evaluation of the angular integration over Im/(p) 
[Eqs. ([T7|), (|19D) and (pOl) ], we find that ImL can be written as a product of Eq. (0) for 
the on-shell two-photon form factor F, times the right hand side of Eq. (|^). 

III. HARMONIC OSCILLATOR WAVE FUNCTION 



We evaluate the quark model form factor R [Eqs. (|l]), (^, (§), (|T2D and (p^T])] using 
a harmonic oscillator parameterization of the qq momentum space wave function that has 



been widely used in nonrelativistic quark model calculations [|T2|,[T3 

$osc(p) = (/?M-'/'exp(^-^). (24) 

Some representative results for the ratio KeR/lmR obtained with this wave function are 
given in Fig. |^. For the sake of illustration, in this figure we take the "physical" mass of 
the hypothetical meson to be equal to the mean energy of the wave packet {mp = 2{Ep)). 

In the limit jS/rriq — ^ 0, the wave packet becomes nonrelativistic, with (Ep) —>■ rriq. 
Bergstrom has made a quark model calculation of R in the extreme nonrelativistic (ENR) 
limit, taking the quark and antiquark to be at rest in the center of mass of the meson, with 
2mg = nip 0. He found: 

Ke(«..,^ia...(i±^)-.„(i±^)4.Sp(-i^)|. (.5, 



Our results agree with Eq. (^) in the hmit l3/mq -^ 0. 

We see that the real part of R is generally small compared to its imaginary part, except 
for small lepton masses mi <C mp, or in the limit of ultrarelativistic binding rUq -C rrip 
(Rei? exhibits a logarithmic divergence as mi or rUq go to zero, which is characteristic of 
mass singularities in loop processes) .0 Since the dispersive part of R appears quadratically 
in the branching ratio, we can expect small corrections to the unitary limit over a wide 
kinematical range. 

To make contact with experimental data on rj —>■ l^l^ , we assign physical values to the 
quark model parameters and meson mass. Various sets of values for the quark masses and 
the wave function Gaussian parameter j3 have been used in the literature [|1^,|1^]. Our results 
for the relative branching fraction Bp are insensitive to variations in these parameters over 
a wide range .0 We use the typical values 



^« = ^TT-d = 330 MeV, 

m, = 550 MeV, (26) 

p = 310 MeV, 

and an rj — 1]' mixing angle of —20° (cf. Ref. [^]). 

We plot the relative branching fraction Bp compared to the unitary bound B^^^ [Eq. 
(I)] in Fig. H, for a range of lepton masses. We take m,p equal to the physical rj mass when 
evaluating Eqs. (|12D and ([T7|). We also include a plot of Bp as obtained by Bergstrom in the 
limit of extreme nonrelativistic binding (cf. Eq. (^) and Ref. 0). The harmonic oscillator 
and ENR results differ by less than 5% over the entire kinematical region in m,i/m,^. 

Using the experimental value for the two-photon branching fraction as input [Q , we find 
B{ri -^ f^'^f^^) = 4.3 X 10~^ and B{r] -^ e^e^) = 6.3 x 10~^ in the harmonic oscillator model. 
The result for the decay to muons is within errors of the recent SATURNE measurement 
Bexpiv ~^ A*^/^ ) = 5.1±0.8 X 10~^ [0, and is only 0.2% larger than the unitary limit. On the 
other hand, the branching fraction to electrons is ~ 3.6 times larger than the corresponding 
unitary limit. We also compute B{7r^ — > e"'"e~) = 1.0 x 10~^ in the oscillator model (about 
twice the unitary limit), although an application of constituent quark models to the pion 
should be viewed with particular caution. 



^We include curves in Fig. |3| to illustrate the results of a harmonic oscillator calculation in the 
ultrarelativistic region; however, the physical region in this model corresponds to f3/mq of order 
unity or smaller. 

^The absolute normalizations of the individual widths T{V -^ 77) and r(P -^ ^u^/i^) depend 
strongly on the choice of model parameters. The widths also depend on the choice of total qq 
energy that is used in the phase space |12,13|. The oscillator wave function results for the two- 



photon widths of the n^, rj and rj' are generally within a factor of a few of the experimental values 

i2i. 



IV. BAG MODEL-INSPIRED RELATIVISTIC WAVE FUNCTION 



Although the mock meson method described in Sec. IIA was developed in the context 
of nonrelativistic quark models |^2|,|13|, we find that the MIT bag model [|l^] can be used to 
motivate a relativistic parameterization of Eq. (]TI1|). We note in this connection that the 
momentum space wave packet distribution $(p) is, in principle, an arbitrary function, and 
need not be localized around momenta that are small compared to the quark mass ||12|| . 

To begin with, we must take account of the fact that bound states in the usual bag 
model are not momentum eigenstates, while the quark and antiquark in Eq. ( pi]| ) for the 
mock meson have equal and opposite three-momentum by construction. Our approach is 
complementary to the formalism introduced by Donoghue and Johnson, who introduced a 
wave packet in order to decompose the bag model wave function into momentum eigenstates 
pO| . Their wave packet is determined by normalizing to a particular matrix element, such 
as /tt. We choose instead to identify the wave packet amplitude $(p) in Eq. (|T0|) directly 
from the Fourier transform of the cavity wave function. 

To begin with, consider the wave function for the ground state of a single quark in a 
cavity of radius R 



^A(r)=Ar 



UJQ 



1/2 



Jl 



j xr\ 



(T-rU 



A 



(r < R), 



(27) 

[xVi?2+^2)l/2_ 

^=R = ig. AT 

= ± along the z 



where for later use we label the "current" quark mass by rrio, and where c^o = 
The momentum eigenvalue x is determined by the boundary condition ipipl 
is a normalization, and U\ is a two-component spinor with polarization A = 
axis. 

The Fourier transform of a Dirac wave function can be found by standard methods ^1 
Since we can decompose the wave function along an arbitrary complete set of states, the 
plane wave spinors in the Fourier transform of Eq. (^) need not have the same ( "current" ) 
quark mass mo as ipxir). We compute a Fourier transform along plane waves of arbitrary 
"effective" mass meg 



^A(r) 



d^ 



P 



1 



(2vr)^/^ ^^[J:^) 



0(p)nA(p)e'P-^ + 0b) E SMp)vy{p)e''''' 

A'=± 



(28) 



where ma(p) and v\{p) are plane wave Dirac spinors with angular momentum A/2 along the 
z axis, and SaaKp) — Ulcr-pUx'- The spinors are normalized to u^u = —v'^v = 2uJes{p), 
where uj'^sip) = p^ + m^^ [||]. 

The Fourier amplitudes and depend on the magnitude of the three momentum 
p= \p\. We find p3| 



<P{p) 



1 



^2tt 



riomo) 



1/2 i?2 /^^g(p) + rneflf 

px \ 



1/2 



X \ Jo{pR -x) - jo{pR + x) + 

p 



2ujcs{p) J 



X- 



' Uesip) +meS 






(29) 



and 



Hp) 



X 



P 



1/2 i?2 / a;^ff(p) + mcff ' 
pa; I 2u;eflf(p) 



1/2 



[JoipR -x) - jo{pR + x)] + 



/uJo-mo\^/'^ 



ujcs{p)+mes \uJo + mo 

X [joipR + x)+ jo{pR -x)- 2jo{pR)jo{x)] \. 



(30) 



The usual mock meson approach to the calculation of hadronic matrix elements assumes 
that the quarks propagate as free particles in intermediate states [0. For example, the 
explicit factor of the quark mass in the matrix elements for V —>■ '~f'~f, l^l^ [Eqs. (^) and 
(p!7|)] is due to the helicity flip along an intermediate quark line. In typical nonrelativistic 
quark model calculations, the quarks have large "constituent" masses. This means, for 
example, that the m, d and s quark components of the rj and rj' flavor wave functions make 
comparable contributions to their matrix elements in these models, in agreement with simple 
phenomenological estimates of the pseudoscalar mixing angle (see, e.g., Ref. 0]). 

In typical MIT bag model calculations however the quarks are assigned current masses, in 
particular, niu^d ~ 0. In order to obtain a sensible phenomenological description of helicity- 
flip amplitudes in a mock meson approach, we must assume that a cavity quark propagates 
with a mass that is different from its "bare" value mo. We therefore assign a constituent 
(or effective) quark mass to intermediate quark lines, and for consistency we use the same 
effective mass in the Fourier transform Eq. (^8]) of the cavity wave function. 

A rigorous calculation of matrix elements in the context of the bag model would use a 
cavity propagator for intermediate quark lines p4|. This was in fact done in a calculation 
of vrO 



77 in Ref. p5 



We note that the cavity quark propagator differs from the free 
propagator by some terms that act, to some extent, like a (momentum-dependent) effective 
mass. A calculation of the leptonic decay using cavity propagators would be quite involved 
however. On the other hand, our extension of the mock meson method to the bag model 
permits a straightforward evaluation of the matrix element incorporating much of the basic 
physics underlying cavity perturbation theory [this is further illustrated by our comments 
below Eq. (|3^) ]. 

Although the actual value of the effective mass for "free" quark propagators is somewhat 
ad hoc, this approach can nevertheless can be used with justification in the calculation of 
the form factor R of interest in this paper [Eq. (|1|)]. This is due to the fact that the ratio 
between the leptonic and two-photon amplitudes is insensitive to variations in the actual 
value of the effective quark mass over a large range. Observe in particular that the same 
factor of the intermediate quark mass due to the helicity flip appears in both matrix elements 
[cf. Eqs. (|T2|) and (|T7p]. On physical grounds, the effective quark mass should be on the 
order of the cavity energy ujo\ we find that Bp changes by only a few percent as TTiefr is varied 
from jujq to ujq. We use meg = ojq in the following calculations. 

At any rate, we do not regard our results as providing precise tests of the predictive power 
of the bag model. Our interest in this paper is to make an estimate of momentum-dependent 
effects in the correction to the unitary limit for V — > l'^l~ in the context of general quark 
models. In that respect, the bag model- inspired wave packet amplitudes provide a very 
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useful comparison with the harmonic oscillator wave function used in Sec. III. In particular, 
(f){p) falls off only as l/p^ at large p, compared to the exponential decay of the oscillator 
wave function. 

In this connection, we note that the mock meson wave packet of Eq. (|l^) is based on 
a single particle description of the quark and the antiquark in the bound state, which is in 
general inadequate to describe localized relativistic states. This problem is shared by the 
single wave packet amplitude introduced by Donoghue and Johnson. A proper connection 
between the localized cavity wave function and momentum space eigenstates requires the 
use of a Bogoliubov transformation [^, which in this case involves a mixing among the 
complete set of cavity eigenstates. 

On the other hand, we find that the negative energy components of the Fourier transform 
actually make a small contribution to the normalization of the ground state cavity mode. 
The cavity wave function Eq. (|27|) is normalized to one in the sphere of radius R, which 
implies 

I d'p {<P\p) + ^\p)) = I. (31) 

We find that the contribution of the negative energy component to this normalization integral 
is at most 7% in the "ultrarelativistic" limit (rrio-R = 0), and falls below 4% by mo-R ~ 1. 
We conclude that a reasonable approximation to the ground state cavity wave function can 
be obtained by truncating the negative energy component of the Fourier transform in Eq. 
(^),P| and we therefore identify $(p) in the mock meson wave packet Eq. ([I0| ) as 

$bag(p) = 0(P) (32) 

[cf. Eqs. ([n]) and (|3)]. 

We note that ReR is very insensitive to the value of the current quark mass mo. This is 
due in large measure to the fact that we use an effective mass nies = 0{ujo) for intermediate 
quark propagators. Thus, even for a massless current quark, the effective mass is nonzero 
{ujqR G [2.04, cxd] for rrioR G [0, cxo]). The logarithmic mass singularity that would be present 
in Re R for a truly massless propagator is thereby avoided. A similar situation occurs in a 
rigorous calculation of the bag model width for tt" -^ 77 using cavity propagators for the 
massless u and d quarks [^ . 

We include our results for the relative branching fraction Bp{ri — > l^l^) for the physical 
7] meson using this bag model-inspired wave function in Fig. §. We use the typical bag 
model parameter values W^ 



0, 

m, = 300 MeV, (33) 

i? = 3.3GeV"^ 

The results for Bp in the relativistic wave function agree with the harmonic oscillator 
calculation (and with the extreme nonrelativistic limit) to within a few percent over wide 



^Note that the normalization of the wave packet cancels in a calculation of the relative branching 
fraction Bp. 
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range in mi/mj^, except for large lepton masses, near threshold. This is evidently related 
to the fact that bag model wave packet has appreciable components at large qq relative 
momentum [cf. <l>bag(p) ~ 1/p^ at large p, compared to $osc(p) ~ exp(— p^)]. 

For the physical cases of rj or vr^ decays to muons or electrons, rrii/mp is small, and the 
bag model-inspired results B{i] -^ /i^yU ) = 4.4 x 10~^, B{i] —>■ e^e^) = 6.1 x 10~^, and 
B{7i^ -^ e+e^) = 1.0 X 10^^ agree with the oscillator model to better than 5%. 

V. SUMMARY 

We have evaluated the electromagnetic box diagram for the leptonic decay of a pseu- 
doscalar quark-antiquark pair with an arbitrary distribution of relative three-momentum. 
Quantitative results were obtained in three different models of the bound state wave func- 



tion (a nonrelativistic harmonic oscillator model |T^JT3[] , a new relativistic momentum space 
wave function that we derived from the MIT bag model, and in the limit of extreme nonrel- 
ativistic binding, analyzed previously by Bergstrom [§). Our results demonstrate that the 
relative branching fraction Bp = B{V — ^ l~^l~)/B{V — >■ 77) is insensitive to the details of 
the quark model wave function. In the case of 77 and 7r° decays, the results obtained with 
the harmonic oscillator and relativistic wave functions agree to within 5%. 

We find that the quark model leads in a natural way to a negligible value for the ratio of 
dispersive to absorptive parts of the electromagnetic amplitude for rj — > fi~^fi~. On the other 
hand, we find substantial deviations from the unitary bound in other kinematical regions, 
such as ?7 — > e"'"e~. 

Using the experimental branching ratio for r] — > 77 as input, these quark models yield 
B{ri —^ P-~^fJ'~) ~ 4.3 X 10^^, within errors of the recent SATURNE measurement of 5.1 ± 
0.8 X 10~^, and B{ri -^ e^e") ~ 6.3 x 10~^. While an application of constituent quark 
models to the pion should be viewed with particular caution, the quark models considered 
here yield B{tt^ — >■ e^e~) ~ 1.0 x 10~^, independent of the details of the model wave 
function to within a few percent. This is to be compared with recent experimental data from 
Fermilab, -Bcxp(7r° -^ e^e~) = 6.9 ± 2.8 x 10~^ [^, and preliminary data from Brookhaven 



-BexplTr" — >■ e~^e~) = 6.0 ± 1.8 x 10~* |]29|. The quark model branching ratios obtained here 
are comparable to the results of a recent analysis in chiral perturbation theory [^, which 
requires the experimental value of the branching ratio for rj -^ IJ'^fJ'" as input. 
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FIGURES 

FIG. 1. Quark model electromagnetic box diagrams for V -^ l^l^ . In a "mock meson" descrip- 
tion (Sec. II), the quark momenta in the center-of-mass are anti-correlated, p,p = (p'^,ibp). The 
lepton momenta are k,k = {k^, itk). 



FIG. 2. Quark model diagrams for V -^ 77. 

FIG. 3. Real part of the form factor R for the relative branching fraction BpiV -^ l~^l~), 
compared to the unitarity result for the imaginary part [cf. Eqs. (|l|) and (P)], in a harmonic 
oscillator model of the bound state wave function [Eq. (p^)]. For this figure, the mass of the 
hypothetical meson is set equal to the mean energy of the qq wave packet. 

FIG. 4. Relative branching fraction Bp for the leptonic decay of the physical r] meson, compared 
to the unitary bound Bp"*. The three curves correspond a harmonic oscillator model [Eqs. (24) 



and (??)], a bag model-inspired relativistic wave function [Eqs. (32) and (??)], and the limit of 
extreme nonrelativistic binding [Eq. ([25|)]. 
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